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The exotic baryon Q + (uudds) is studied with microscopic calculations in a quark model by 
using a method of antisymmetrized molecular dynamics. We predict that three narrow states, 
,F = 1/2+ (I = 0), J* = 3/2+ (J = 0), and J w = 3/2" (I = 1) nearly degenerate with the lowest 
l/2~ state in the uudds system. We discuss KN decay widths and estimate them to be V < 7 for 
the J n = {l/2 + , 3/2 + }, and F < 1 MeV for the J w = 3/2~ state. In contrast to these narrow states, 
the l/2~ states should be much broader. We assign the observed 9+ as the ,F = {1/2+ , 3/2+}. 
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The exotic baryon + has recently been reported by several experimental groups [1-9]. Since the quantum numbers 
determined from its decay modes indicate that the minimal quark content is uudds, these induced experimental and 
theoretical studies of multiquark hadrons. However it should be kept in mind that the + has not been well established 
yet because of the low statistics and experimental reports [10-12] for no evidence of the + . 

The prediction of a J x = l/2 + state of uudds by a chiral soliton model [13] motivated the experiments of the first 
observation of + [1]. Their prediction of even parity is unnatural in the naive quark model, because the lowest q q 
state is expected to be spatially symmetric and have odd parity due to the odd intrinsic parity of the anti-quark. 
Theoretical studies were done to describe + by many groups [14-21], some of which predicted the opposite parity, 
= 1/2 - [18-20]. The problem of spin and parity of + is not only open but also essential to understand the 
dynamics of pentaquark systems. To solve this problem, it is crucial to calculate five-quark system relying on less a 
■ priori assumptions such as the existence of quark clusters or the spin parity. 

In this paper we would like to clarify the mechanism of the existence of the pentaquark baryon and predict possible 
narrow + states. We try to extract a simple picture for the pentaquark baryon with its energy, width, spin, parity 
and also its shape from explicit 5-body calculation. In order to achieve this goal, we study the pentaquark with a 
Qh, flux-tube model [24,25] based on strong coupling QCD, by using a method of antisymmetrized molecular dynamics 
(AMD) [22,23]. In the flux-tube model, the interaction energy of quarks and anti-quarks is given by the energy of 
\ the string-like color-electric flux, which is proportional to the minimal length of the flux-tube connecting quarks and 
anti-quarks at long distances supplemented by perturbative one-gluon-exchange (OGE) interaction at short distances. 
For the q 4 q system the flux-tube configuration has an exotic topology, Fig. 1(c), in addition to an ordinary meson- 
baryon topology, Fig. 1(d), and the transition between different topologies takes place only in higher order of the 
strong coupling expansion. Therefore, it seems quite natural that the flux-tube model accommodates the pentaquark 
baryon. In 1991, Carlson and Pandharipande studied exotic hadrons in the flux-tube model [26]. They calculated for 
only a few q A q states with very limited quantum numbers and concluded that pentaquark baryons are absent. We 
apply the AMD method to the flux-tube model. The AMD is a variational method to solve a finite many-fermion 
system. This method is powerful for the study of nuclear structure. One of the advantages of this method is that the 
spatial and spin degrees of freedom for all particles are independently treated. This method can successfully describe 
various types of structure such as shell-model-like structure and clustering (correlated nucleons) in nuclear physics. In 
the application of this method to a quark model, we take the dominant terms of OGE potential and string potential 
due to the gluon flux tube. Different flux-tube configurations are assumed to be decoupled. Since we are interested 
in the narrow states, we only adopt the confined configuration given by Fig. 1(c). We calculate all the possible spin 
parity states of uudds system, and predict low-lying states. By analysing the wave function, we discuss the properties 
of + and estimate the decay widths of these states with a method of reduced width amplitudes. 

This paper is organized as follows. We explain the formulation of the present framework in the next section, and 
show the results in III. In IV, we discuss the structure of low-lying states and their widths. Finally, we give a summary 
in V. 
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II. FORMULATION 



In the present calculation, the quarks are treated as non-relativistic spin-i Fermions. We use a Hamiltonian as 
follows, 

H = H + Hi + Hf, (1) 

where H is the kinetic energy of the quarks, Hj represents the short-range OGE interaction between the quarks and 
Hf is the energy of the flux tubes. For simplicity, we take into account the mass difference between the ud quarks 
and the s quark, only in the mass term of Hq but not in the kinetic energy term. Then, Hq is represented as follows; 



ffo = £m i + £^--T , (2) 



where N q is the total number of quarks and is the mass of i-th quark, which is m q for a u or d quark and m s for 
a s quark. To denotes the kinetic energy of the center-of-mass motion. 

Hi represents the short-range OGE interaction between quarks and consists of the Coulomb and the color-magnetic 
terms, 



Hi = a c J2 F i p j 
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Here, a c is the quark-gluon coupling constant, and FiFj is defined by X^ a =i 8 F^F", where Ff is the generator of 
color SU(3), \\f for quarks and — \(\f)* for anti-quarks. The usual S(rij) function in the spin-spin interaction is 
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as in Ref. [26]. Of course, the full OGE interaction 



basic properties of the pentaquark, we do not include these minor contributions. 

In the flux-tube quark model [24], the confining string potential is written as Hf = aLf — M°, where a is the 
string tension, Lf is the minimum length of the flux tubes, and M° is the zero-point string energy. M° depends 
on the topology of the flux tubes and is necessary to fit the qq, q 3 and q 4 q potential obtained from lattice QCD or 
phenomenology. In the present calculation, we adjust the M° to fit the absolute masses for each of three-quark and 
pentaquark. 

For the meson and 3g-baryon systems, the flux-tube configurations are the linear line and the Y-tyj>e configuration 
with three bonds and one junction as shown in Fig. 2(a) and (b), respectively. The string potential given by the F-type 
flux tube in a 3g-baryon system is supported by Lattice QCD [27]. For the pentaquark system, the different types of 
flux-tube configurations appear as shown in Fig. l.(e),(f), and (d), which correspond to the states, |$( e )) = I [ud] [ud] s) , 
[$(/)) = |[uu][cZd]s), and | = \(qqq)i(qq)i) , respectively, ([qq] is defined by color anti-triplet of qq.) The flux- 
tube configuration (e) or (f) have seven bonds and three junctions, while the configuration (d) has four bonds and 
one junction. In principle, besides these color configurations ([qq] [qq]q and (qqq)i(qq)i), other color configurations 
are possible in totally color-singlet q A q systems by incorporating a color-symmetric (qq)% pair as in Refs. [16,21]. 
However, since such a string from the (qq)e is energetically excited and is unfavored in the strong coupling limit of 
gauge theories as shown in Ref. [28]. Therefore, we consider only color-3 flux tubes as the elementary tubes. In 
fact, the string tension for the color-6 string in the strong coupling limit is 5/2 times larger than that for the color-3 
string from the expectation value of the Casimir operator. The string potentials given by the tube lengths of the 
configuration Fig. 1(c) is supported by Lattice QCD calculations [29]. 

In the present calculation of the energy, we neglect the transition among |$( e )}, |$(/)) an d |3>(d)) because they have 
different flux-tube configurations. It is reasonable in the first order approximation, as mentioned before. In each tube 
configuration, the minimum length Lf is given by a sum of the lengths (Ri) of bonds Lf = R\ + ■ ■ ■ + (k is the 
number of the bonds. See Fig. 2). Here we define Lij to be the length of the path between i-th (anti)quark and j-th 
(anti)quark along the flux tubes. For example, in case of the [gqllgq]? state shown in Fig. 2(c), the path lengths are 
given by the bond lengths Ri as L\ 2 = R\ + i?2, £-13 = Ri + Re + Ri + R3, L x \ = Ri + i? 6 + -R5, etc. Then we 
can rewrite Lf in the expectation values of the string potential ($|ii/|<I>) with respect to a meson system($ g q), a 
three-quark system^^), and the pentaquark states <f>( e ), $(/), as follows: 

L f = L 12 in <* M -|If/l*«->. (4) 

L f = ^(L 12 +L 2 3 + L 31 ) in (<$> q3 [Hf\<$> q 3), (5) 
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FIG. 1. Flux-tube configurations for confined states of qq (a), q" 
(e) and (f) represent the flux tubes in the color configurations, [uc 



(b), q 4 q (c), and disconnected flux-tube of q 4 q (d). Figures 
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In the practical calculation, we approximate the minimum length of the flux tubes Lf by a linear combination of 
two-body distances m between the i-th (anti)quark the j-th (anti)quark as, 



L f w n 2 in ($^1^1$^), 

~ t^ 7 " 12 + r 23 + r 3 i) in (^^l^l*^), 

£/ ~ ^(ri2 +r 34 ) + ^(ri3 + r 14 + r 23 + r 24 ) + ^(r Tl + r l2 + r l3 + r l4 ) in (^ (eJ) \H f \^ {eJ) ), 
1 



(8) 
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Lf « -(ri2 +r 23 +r 3l ) + r l4 in (S^iff |$ (d) ). 

It is clear that the above equations are obtained by approximating the path length L^ with the distance as L^ w 
for all gg and qq pairs. In the meson system, it is clear that Eq.8 gives the exact Lf value. The approximation, Eq.9, 
for 3g-baryons is used in Ref. [24] and has been proved to be a good approximation. We note that the confinement 
is reasonably realized by the approximation in Eq.10 for &( e ,f) as follows. The flux-tube configuration (e)(or (f)) 
consists of seven bonds and three junctions. In the limit that the length(i?i) of any i-th bond becomes much larger 
than other bonds, the string potential (Hf) approximated by Eq.10 behaves as a linear potential oR. It means that 
all the quarks and anti-quarks are bound by the linear potential with the tension a. In that sense, the approximation 
in Eq.10 for the connected flux-tube configurations is regarded as a natural extension of the approximation(Eq.9) for 
3<7-baryons. It is convenient to introduce an operator O = — j<?J2i<j FiFj^ij — M°. One can easily prove that the 
above approximations, 8,9,10,11, are equivalent to ss within each of the flux-tube configurations 

because the proper factors arise from FiFj depending on the color configurations of the corresponding qq (or qq) pairs. 

In order to see the accuracy of the approximations Eqs.9 and 10, we calculate the ratio of the approximated length 
Lapp to the exact Lf in a simple quark distribution with Gaussian form which imitates the model wave function of 
the present calculation. Figure 3 shows the ratio L app /Lf in a q 3 system and a [qq] [qq] q system. The quark positions 
Yi are randomly chosen in Gaussian deviates with the probability p = exp(— rf /b 2 ), and (Lf, L app /Lf) values for 1000 
samplings are plotted. We use the same size parameter b as that of the single-particle Gaussian wave function in 
the present model explained later. Comparing Figs. 3(a) with 3(b), we found that the L app /Lf ratio for the [qq] [qq]q 
system is about 10% smaller than that for the q 3 system. Since the zero-point energy Mq in the string potential is 
adjusted in each of the q 3 and the [qq] [qq]q, this underestimation should relate only to the relative energy of the string 
potential in each system, and may give a minor effect on the level structure of the pentaquark. 

We solve the eigenstates of the Hamiltonian with a variational method in the AMD model space [22,23]. We take 
a base AMD wave function in a quark model as follows. 



*(Z) = (1±P)A [<f>z 1 <Pz 2 ---4>z Nq X\ , 
. = (^) 3/4 exp[-^(r-V26Z, t ) 2 + iZf], 



(12) 
(13) 
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FIG. 2. Flux-tube topologies for the qq (a), q 3 (b), [9192] [g3<?4]<7i (c), and disconnected flux tubes (d) for the (qqq)i(qq)i. 
The flux-tube topologies are described by the bonds with the lengths Rk and the junctions ju- 
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FIG. 3. The ratio L apP /Lf for the approximated tube length L app and the exact tube length Lf in a g 3 system and a 
[qq][qq]q system. The quark positions are randomly chosen in Gaussian deviates with the probability p — exp(— r 2 /b 2 ), and 
(Lf, L app /Lf) values for 1000 samplings are plotted. 
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where 1 ± P is the parity projection operator, A is the anti-symmetrization operator, and the spatial part </>Zi of the 
i-th single-particle wave function given by a Gaussian whose center is located at Z^ in the phase space. X is the spin- 
isospin-color function. For example, in case of the proton, X is given as X = (| tlt)s ~~ I TT4)s) ® \uud) ® e a bc\dbc)c- 
Here, \m)s{m =T,I) is the intrinsic-spin function and \a)c{a = 1,2,3) expresses the color function. Thus, the wave 
function of the N q quark system is described by the complex variational parameters, Z = {Zi, Z 2 , • • • , Zjv,}- By using 
the frictional cooling method [22] the energy variation is performed with respect to Z. 
For the pentaquark system (uudds), 

X = c mim2m3mim5 \m 1 m 2 m 3 rn i m 5 )s ® {\ududs} or \uudds}} <g> e 

abg^ceh^ghf \abcef) c , (14) 



where \ududs) and \uudds) correspond to the configurations [ud][uc?]s and [uu][c?d]s in Fig.l, respectively. Since 
we are interested in the confined states, we do not use the meson-baryon states, (qqq)i(qq)i- This assumption of 
decoupling of the reducible and irreducible configurations of the flux tubes can be regarded as a kind of bound-state 
approximation. The decoupling of the different flux-tube configurations can be characterized by the suppression 
factor e from the transition of the gluon field in the non-diagonal matrix elements e($i |0|<f> 2 )- In a simple flux-tube 
model, e is roughly estimated by the area A s swept by the tubes when moving from one configuration into the other 
configuration as e ~ exp(— aA s ). We make an estimation of the expectation value of exp(— aA s ) by assuming a simple 
quark distribution with Gaussian form which imitates the model wave function in the same way as the evaluation of 
the L app /Lf. The suppression factor e among the configurations [u<i][u<i]s, [mt][<ic/]s, and {qqq)i{qq)i is estimated to 
be e 2 < 1/10 within the present model space. Therefore, we consider that the present assumption of the complete 
decoupling e = in the energy variation is acceptable in first order calculations. 

The coefficients c mim2m3mims for the spin function are determined by diagonalization of Hamiltonian and norm 
matrices. After the energy variation with respect to the {Z} and c mim2m3m4rriB , the intrinsic-spin and parity S n 
eigen wave function $(Z) for the lowest state is obtained for each S v . In the AMD wave function, the spatial wave 
function is given by multi-center Gaussians. When the Gaussian centers are located in some groups, the wave function 
describe the multi-center cluster structure and is equivalent to the Brink model wave function(a cluster model often 
used in nuclear structure study) [30,31]. On the other hand, because of the antisymmetrization, it can also represent 
shell-model wave functions when all the Gaussian centers are located near the center of the system [30,31]. In nuclear 
structure study, it has been already proved that the AMD is one of powerful tools due to the flexibility of the wave 
function [23]. In general, the relative motions in the AMD are given by such Gaussian forms as exp[— z/(x — R) 2 ] 
where x is a Jacobi coordinate and R is given by a linear combination of the Gaussian centers {Zi, Z 2 , • • • , Zjv,}- 
Here we explain the details of the relative motion in a simple case of a two-body cluster structure in a N q = 5 system. 
If the Gaussian centers are located in two groups as Zi = Z 2 = Z3 = Qi/\/2b and Z4 = Z5 = Q 2 /\/2fr, and if each 
group does not contain identical particles, the wave function expresses the two-body cluster state, where each cluster 
is the harmonic oscillator Os-orbital state, (0s) 2 ' 3 , with zero orbital-angular momentum. The inter-cluster motion X 
is given as X(x, R, v') = exp[— z/(x — R) 2 ], where v' = g 3 ^, R = Q 2 Qi and x is the relative coordinate between 
the clusters. In the partial wave expansion of the inter-cluster motion X, 



where i L is the modified spherical Bessel function, it is found that the wave function contains higher orbital-angular 
momentum L components in general. However, in case of v' ' R 2 < 0(1), the wave function is dominated by the lowest 
L component since the L components rapidly decrease with the increase of L. As a result, the even-parity S 1T=+ 3 q 
and odd-parity q A q states are almost the L = eigen states, while the odd-parity 3 9 and even-parity q A q states are 
nearly the L — 1 eigen states. (The q 4 q contains an odd intrinsic parity of the q in addition to the parity of the spatial 
part.) Therefore, we do not perform the explicit L-projection in present calculation for simplicity. We have actually 
checked that the obtained wave functions are almost the L-eigen(L = or 1) states and higher L components are 
minor in most of the q 3 and q q states. 

In the present wave function we do not explicitly perform the isospin projection, however, the wave functions 
obtained by energy variation are found to be approximately isospin-eigen states in most of the low-lying states of the 
q 3 and q 4 q due to the color-spin symmetry. 

In the numerical calculation, the linear and Coulomb potentials are approximated by seven-range Gaussians. We 
use the following parameters, 




X(x, R, v') = cxp[-i/(x - R) 2 ], 
= J2l ^ L {2v'Rx)e-»'^+ R2 ) Em Y lm (x)Y* m (R), 



(15) 



a c = 1.05, 
A = 0.13 fm, 
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TABLE I. Calculated masses (GeV) of the q 3 systems. The expectation values of the kinetic, string, Coulomb and 
color-magnetic terms are also listed. 
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TABLE II. Calculated masses(GeV) of the uudds system. M° 4 _=2385 MeV is used to adjust the energy of the lowest state to 
the observed mass. The expectation values of the kinetic, string, Coulomb, color-magnetic terms, and that of the color-magnetic 
term in qq pairs are listed. In addition to the lowest l/2~ state with the [uu][dd]s configuration, we also show the results of 
the l/2~ state with [iid][iid]s configuration, which lies in the low-energy region. 
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m q = 0.313 GeV, 
m s = 0.513 GeV, 

a = 0.853 GeV/fm. (16) 

Here, the quark-gluon coupling constant a c is chosen so as to fit the N and A mass difference. The string tension a 
is adopted to adjust the excitation energy of AT* (1520). The size parameter b is chosen to be 0.5 fm. 

III. RESULTS 

In tabled, we display the calculated energy of q 3 states with S n = 1/2+ (N), = 3/2+ (A), S n = 1/2~(N*). The 
zero-point energy M° of the string potential is chosen to be M° 3 = 972 MeV to fit the masses of q 3 systems, N, N* 
and A. The calculated masses for A with S n = l/2~ and l/2 + correspond to the experimental data of A(1115) and 
A* (1670). The contributions of the kinetic and each potential terms are consistent with the results of the Ref. [26]. 
We checked that the obtained states are almost eigen states of the angular momentum L and the L projection gives 
only minor effects on the energy. 

Now, we apply the AMD method to the uudds system. For each spin parity, we calculate energies of the [wd][ud]s 
and [uu][c?d]s states and adopt the lower one. In table. II, the calculated results are shown. We adjust the zero-point 
energy of the string potential M as M® 4 - = 2385 MeV to fit the absolute mass of the recently observed 9 + . This 
M° 4 - for pentaquark system is chosen independently of M° 3 for 3g-baryon. If Af° 4 - = §M° 3 is assumed as Ref. [26], 
the calculated mass of the pentaquark is around 2.2 GeV, which is consistent with the result of Ref. [26]. 

The most striking point in the results is that the S 11 = 3/2~ and = 1/2+ states nearly degenerate with the 
S n = 1/2- states. The S* = 1/2+ correspond to J* = {l/2+,3/2+} (S = 1/2, L = 1), and the S n = 3/2" is 
J v = 3/2~ (S = 3/2, L = 0). The lowest S w = l/2-(J 7T = l/2~, L = 0) state appears just below the = 3/2" and 
the second S n = l/2"(J 7r = l/2~, L = 0) state is at the same energy as the S 7 ' = l/2+(J 7r = l/2+,3/2+,L = 1) 
states. However these J v = l/2~ states, as we discuss later, are expected to be much broader than other states. 
The J 71 = 1/2+ and 3/2+ exactly degenerate in the present Hamiltonian which does not contain the spin-orbit force. 
Other spin-parity states are much higher than these low-lying states. 
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IV. DISCUSSION 



In this section, we analyze the structure of the obtained low-lying states of the uudds system, and discuss the level 
structure and the width for KN decays. 

A. Structure of low-lying states 

We analyze the spin structure of these states, and found that the J 77 = {l/2+,3/2+} states consist of two spin-zero 
[ud]-pairs, while the J 77 = 3/2" contains of a spin-zero [wcij-pair and a spin-one [w<i]-pair. Here we call the color 
anti-triplet qq pair with the same spatial single-particle wave functions as a [qq]-pair and note a spin S [qqj-pair 
as [qq] s ■ Since the [ud]°-pair has the isospin 1 = and the [wd] 1 -pair has the isospin 1=1 because of the color 
asymmetry, the J 77 = 3/2" state is isovector while the lowest even-parity J 77 = {l/2+,3/2+} states are isoscalar. 
The J 77 = 1/2+ state corresponds to the 6 + (1530) in the flavor 10-plet predicted by Diakonov et al. [13]. It is 
surprising that the odd-parity state, J 77 = 3/2" has the isospin 1=1, which means that this state is a member of 
the flavor 27-plet and belong to a new family of 9 baryon. We denote the J 71 ' — {1/2+, 3/2+}, / = states by 0~| , 
and the J 77 = 3/2", I = 1 state by Of. The mass difference E(G^) - E(Of) is about 30 MeV. In the energy region 
compatible to the J 71 " = {l/2 + ,3/2 + } and J 71 " = 3/2~ states, there appear two J 71 " = {1/2"} states. The lowest one 
is the [liujfddjs state with [mm] 1 and [dd] 1 pairs, while the higher one is the [mc^uc 7 ] 7 ? with [ud]° and [ud] 1 pairs. The 
former is the isospin symmetric state and is dominated by / = component. The latter is isovector and is regarded 
as the spin ^-partner of the J 77 = 3/2" state. The J 71 ' = 1/2" state is the lowest in the uudds system. We, however, 
consider this state not to be the observed + because its width should be broad as discussed later. 

Although it is naively expected that unnatural spin parity states are much higher than the natural spin-parity 1/2" 
state, the present results show the abnormal level structure of the (ududs) system, where the high spin J 77 = 3/2~ 
state and the unnatural parity J 71 " = {l/2 + ,3/2 + } states nearly degenerate just above the J 77 = 1/2" state. By 
analysing the details of these states, the abnormal level structure can be easily understood with a simple picture 
as follows. As shown in table. II, the J 77 = {1/2+, 3/2 + }(L = 1) states have larger kinetic and string energies than 
the J 77 = 3/2~(L = 0) and J 77 = 1/2" (L = 0) states, while the former states gain the color-magnetic interaction. 
It indicates that the degeneracy of parity-odd states and parity-even states is realized by the balance of the loss of 
the kinetic and string energies and the gain of the color-magnetic interaction. In the J 77 = {l/2 + ,3/2 + } and the 
J 77 = 1/2", 3/2" states, the competition of the energy loss and gain can be understood by Pauli principle from the 
point of view of the [ggj-pair structure as follows. As already mentioned by Jaffe and Wilczck [14], the relative motion 
between two [<7q]°-pairs must have the odd parity (L = 1) because the L = is forbidden between the two identical 
[ggj-pairs due to the color antisymmetry. In the J 77 = 3/2" state and the second J 77 = 1/2" state, one of [ud]°-pairs is 
broken to be a [wd] 1 -pair and the L = is allowed because two diquarks are not identical. The L = is energetically 
favored in the kinetic and string terms, and the energy gain cancels the color- magnetic energy loss of a [uc 7 ] 1 -pair. 
Also in the lowest J 77 = 1/2" state, the competition of energy loss and gain is similar as each contribution of the 
kinetic, string and potential energies in the lowest J 77 = 1/2" state is almost the same as those in the J 77 = 3/2" and 
the second J 77 = 1/2" (table II). It means that the gain of the kinetic energy of the L = state compete with the 
color-magnetic energy loss in the lowest J 77 = 1/2" as well as the J 77 = 3/2" and the second J 77 = 1/2". 

We should stress that the existence of two spin-zero ud-diquarks in the J 77 = {l/2+,3/2+} states predicted by 
Jaffe and Wilczek [14] is actually confirmed in the present calculations without a priori assumptions for the spin and 
spatial configurations. In fact, the component with two spin-zero [wd]-pairs is 97% in the present J 77 = {l/2 + ,3/2 + } 
state. In Fig. 4, we show the quark and anti-quark density distribution in the J 77 = {l/2 + ,3/2 + } states and display 
the centers of Gaussians for the single-particle wave functions. In the intrinsic wave function, Gaussian centers for 
two [u(i] -pairs are located far from each other with the distance about 0.6 fm. It indicates the spatially developed 
diquark-cluster structure, which means the spatial and spin correlations in each [udj-pair. It is found that the center 
of the s stays at the same point of that of one [ud]° 7 as Zi = Z 2 = 7-75^12 and Z 3 = Z 4 = Z 5 = — t~t^Qi2 where 
{Zi, Z2, • • • , Z5} are the Gaussian centers in Eq.12 and IQ12I ~ 0.6 fm. As a result, we found the spatial development 
of ud-uds clustering and a parity-asymmetric shape in the intrinsic state before parity projection(Fig.4). As explained 
in II, the wave function is equivalent to the [ud} -[ud\°s cluster wave function in Brink model [30] with L = 1 relative 
motion. After the parity projection, the s is exchanged between two diquarks. In contrast to the spatially developed 
cluster structure in the even-parity state, the odd-parity states J 77 = 1/2", 3/2" are almost the spatially symmetric 
(Os) 5 states with spherical shapes. 

As mentioned before, the degeneracy of the even-parity states and the odd-parity states originates in the balance 
of the L = 1 excitation energy and the energy gain of the color-magnetic interaction. Here we consider the L = 1 
excitation energy AE(L = 1) as the total energy lose in the kinetic, string and Coulomb terms. It is important that 
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FIG. 4. The q and q density distribution in the J" = 1/2+ , 3/2+ (5 1 = 1/2, L = 1) states of the uudds system. The u density 
(a), s density (b), and total quark-antiquark density (c) of the intrinsic state before parity projection are shown. The schematic 
figure of the corresponding flux-tube configuration is illustrated in (d). Open squares in (a) and (b), indicate the positions of 
Gaussian centers Re[\/2&ZJ for the i-th single-particle wave functions. 



AE(L = 1) <~ 0.3 GcV in the pentaquark is much smaller than AE(L = 1) <~ 0.5 GeV in the nucleon system. The 
reason for the relatively small AE{L = 1) in the pentaquark can be easily understood by the ud-uds cluster structure. 
In the two-body cluster state with the L = 1 relative motion, the AE(L = 1) is roughly estimated by the reduced 
mass u = AiA 2 /(Ai + A 2 ) of two clusters, as is given as AE(L = 1) oc - (Ai and A 2 are the masses of the clusters). 
In the nucleon, u = |m g is obtained from the ud-u cluster structure in the J" = 1/2~(L = 1) state, while u <~ |m g 
for the pentaquark system is found in the ud-uds clustering. The reduced mass in the pentaquark is 9/5 times larger 
than that in the nucleon system, therefore, AE(L = 1) should be smaller in the pentaquark than in the nucleon by 
the factor 5/9. This factor is consistent with the present AE(L = 1) values. 

We give a comment on the LS'-splitting between J v = 1/2+ and 3/2 + (S l = 1/2, L = 1). In the present calculation, 
where the spin-obit force is omitted, the J w = 1/2+ and 3/2 + states exactly degenerate. Even if we introduce the 
spin-orbit force into the Hamiltonian, the LS'-splitting should not be large in this diquark structure because the effect 
of the spin-orbit force from the spin-zero diquarks is very weak as discussed in Ref. [33] . 

We remark that the [ud]°-[ud]°s cluster structure in the present result is different from the diquark-triquark structure 
proposed by Karliner and Lipkin [16] because the wcis-triquark in Ref. [16] is the (us)g =1 s with the color-symmetric 
spin-one ue?-diquark. In the (ws)g =1 s-triquark, the s quark should be tightly bound in the triquark due to the strong 
color-magnetic interaction between (us)f =1 and s. On the other hand, in the present [ud] °s-cluster, the s feels no 
strong color-magnetic interaction and is bound more weakly than in the (ws)f =1 s-triquark. The color-6 flux tubes 
are not taken into account in the present framework since they are excited. However, the (ws)g =1 s-triquark might be 
possible if the short-range correlation in the triquark make the flux-tube short enough to be excited into the color-6 
flux-tube. 



B. Width for KN decays 

In the 6 + — > KN decays, it is important that the allowed decay mode in the Of (J 71 " = 3/2~) is D wave, which 
should make the Of state narrower than the Oq (J v = 1/2+, 3/2+) because of higher centrifugal barrier. We estimate 
the ifiV-decay widths of these states by using a method of reduced width amplitudes [31,32]. This method has been 
applied for the study of cv-decay width in the nuclear physics within bound state approximations. In this method, the 
decay width T is estimated by the penetrability o) of the barrier and the reduced width 7 2 (a) as a function of 

the threshold energy E t h and the channel radius a, 



P L {k,a) - 



r = 2P L (fc,a) 7 2 (a), 
ka 



j 2 L (ka) + n\{ka) 
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7 2 («) = |^ /oc (a), (17) 



where n is the reduced mass, k is the wave number k = y 2/iE t h/?i 2 , and j'i,(ti£,) is the regular (irregular) spherical 
Bessel function. Sf ac (a) is the probability of decaying particle at the channel radius a. We define Y a L (a,E t h) = 
^JT j'i (ka)+n' A (ka) ' tncn ; the decay width can be rewritten in a simple form as T = x Sf ac . In the following 
discussion, we choose the channel radius a = 1 fm and E th = 100 MeV. Since the transitions between the different 
flux-tube configurations, a confined state [«.d][u(i]s and a decaying state (udd) i (us) i , are of higher order, the Sf ac 
should be small in general when the suppression by the flux-tube transition is taken into account. Here, we evaluate 
the maximum values of the widths for the J 7 ' = 1/2+ , 3/2+ states with the method of the reduced width amplitudes, 
by using meson-baryon probability considering only the simple overlap for the quark wave functions. 

In case of even-parity J 7 ' ~ l/2 + ,3/2 + states, the KN decay modes are the P-wave, which gives « 100 

MeV fm. By assuming (Os) 2 and (0s) 3 harmonic-oscillator wave functions for K° and p, we calculate the overlap 
between the obtained pentaquark wave function and the K°p state. As explained in the previous subsection, the J 71 " = 
l/2+,3/2 + states have the ud-uds cluster structure where five Gaussian centers are written as Zi = Z 2 = ■g^;Qi2 

and Z 3 = Z4 = Z 5 = — ^j^Qi2- We assume a simple K°p wave function as follows, 

$if o p = (1 + P)A [<£ Zl 0z 2 ■ ■ ■ H Nq x] , (18) 

<fe = (^) 3/4 cxp [-^(r - V2bZ t r + iZf] , (19) 

where the Za are chosen as Zi = Z 2 = Z 3 = a^ 2 ^Qi 2 /|Qi2|, Z 4 = Z 5 = — a^7^Qi 2 /|Qi 2 |, and the spin-isospin-color 
wave function is taken to be 

X= c mim2m3m4ms \m 1 m 2 m 3 m4m 5 } s ®\ududs} ®e abc 5 ef \abcef)c- (20) 

mi,m2,m3,m4,m5 

The same size parameter b as that of the pentaquark is used. The coefficients c mim2m3m4TOs for the spin function 
are taken to express the J 71 = 1/2+ proton and the pseudoscalar K° meson. The probability Sf ac is evaluated by 
the overlap with the obtained J 71 " = l/2+,3/2+ wave function, Sf ac = |($^-o p |$(Z))| 2 . (The <& K o p and $(Z) arc 
normalized.) The probability Sf ac = 0.034 fm -1 is evaluated by the overlap. Roughly speaking, the main factors in 
this meson-baryon probability are the factor | from the color configuration, the factor j from the intrinsic spin part, 
and the other factor which arises from the spatial overlap. By using the probability Sf ac = 0.034, the K p partial 
decay width is evaluated as T < 3.4 MeV. The K + n decay width is the same as that of the K°p decay, and the total 
width of the J 71 " = l/2 + , 3/2+ states is estimated to be T < 7 MeV. This is consistent with the discussion in Ref. [34]. 

It is interesting that the KN decay width of the @f( J T = 3/2~) is strongly suppressed by the D-wave centrifugal 
barrier, because lower spin (5-wave and P-wave) decays are forbidden due to the conservation of spin and parity. 
Consequently, T° =2 is w 30 MeV fm, which is much smaller than that for the P-wave case. Moreover, the 0f (J 17 = 
3/2~) is the state with S 71 = 3/2~ and L = 0, which has no overlap with the KN(S 7T = l/2~ and L = 2) states in the 
present calculation because the spin-orbit or tensor forces are omitted. If we introduce the spin-orbit or tensor forces, 
the D-state(S 7T = l/2~ and L = 2) will be slightly mixed into the (J 71 = 3/2 _ ). However, the mixing component 
should be small because of the dominant central force in the potential. In other words, the KN probability(S'/ oc ) 
in the ^(J 71 = 3/2~) state is expected to be rather suppressed than that in the Q^iJ 71 = l/2+,3/2+) states. 
Considering the suppression effects in both terms T° and <S/ ac , the J n — 3/2~ state should be extremely narrow. 
If we assume the Sf ac in the J 7 ' = 3/2~ to be half of that in the J 71 = l/2 + ,3/2 + states, the KN decay width is 
estimated to be T < 1 MeV. 

Contrary to the narrow width of the J 71 = 3/2~ state, the J 71 = l/2~ state should be much broader than other 
states because S'-wave(i = 0) decay is allowed and therefore the centrifugal barrier is absent. We cannot evaluate the 
width of the J 71 = l/2~ states with the present method, since the method of the reduced width amplitudes works only 
when there exist barriers in the decaying channels. If we adopt the theoretical width T = 1.1 GcV for the J" = l/2~ 
states in [34] and the suppression factor e 2 < 1/10 due to the string transition, the width is evaluated to be T <~ 100 
MeV, which is still too large to describe that of the observed + . We consider that the J 71 = l/2~ states may melt 
away due to the coupling with KN continuum states with no centrifugal barrier. 

Also in other spin-parity states, the coupling with the KN continuum states is important for more quantitative 
discussions on the widths. We should point out that, in introducing the meson-baryon coupling, one should not 
treat only the quark degrees of freedom but take into account the suppression due to the rearrangement of flux-tube 
topologies between the meson-baryon states and the confined states. 
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V. SUMMARY 



We proposed a quark model in the framework of the AMD method, and applied it to the uudds system. The level 
structure of the uudds system and the properties of the low-lying states were studied within the model space of the 
[iq] [<1Q]q configuration, where all the (anti)quarks are connected by the color-3 flux tubes. We predicted that the 
narrow J 71 = {l/2 + ,3/2 + }(9 ) and J 77 = 3/2" (6i) states nearly degenerate with the J 77 = 1/2" states. The widths 
of the J 71 = {1/2+, 3/2+} states and the 3/2" state are estimated to be T < 7 MeV and T < 1 MeV, respectively. On 
the other hand, the J 7 ' = l/2~ states should be broad, and we consider that they may melt away due to the coupling 
with KN continuum states with no centrifugal barrier. Two spin-zero diquarks are found in the {l/2 + , 3/2+} states, 
which confirms Jaffc-Wilczek picture. We comment that the formation of two spin-zero diquarks does not always 
occur in J" = {1/2+, 3/2+} pentaquarks. For example, in case of the ddssu system, the diquark structure disappears. 
Instead, a dss-du clustering appears in the J OT = {l/2 + , 3/2 + } [ds][c?s]w states because the color magnetic interaction 
is weaker for ds pairs than for du pairs in OGE potential. In other words, the diquark structure is formed in such 
a certain pentaquark as the Oq due to the strong color-magnetic attraction between ud quarks. The degeneracy of 
the J 77 = 1/2+ 3/2+ 1/2+ and 3/2+ states is realized by the balance of the kinetic and string energies and the 
color-magnetic interaction. The origin of the novel level structure is due to the color structure in the confined five 
quark system bound by the connected flux-tubes. 

The J 77 — {l/2 + , 3/2 + }(9j =0 ) states in the present results may be assigned to the experimentally observed © + , 
while J 77 = 3/2 _ (6/ = i) is not observed yet. One should pay attention to the properties of these states, because the 
production rates depend on their spin, parity and widths. The existence of many narrow states, J 77 = 1/2+, 3/2 + , 
and 3/2+ for the Oq and may help to explain the inconsistent mass positions of the 9 + among the different 
experiments. Especially, the double peaks of the J 77 = 1/2+ and 3/2 + states in the ©J are expected. In the 0~]~ 
peaks observed in the invariant mass or missing mass spectra, it is difficult to find the possible double peaks because 
the statistics and the resolutions are not enough [19]. The analysis of the NK scattering [35] provided the upper 
limit T < 1 MeV for the widths of each peaks. Considering the suppression factor due to the gluon transitions, the 
possibility of the double peaks (J 71 " = 1/2+ and 3/2+) suggested in the present works has not been excluded yet. We 
should comment that another explanation for the inconsistency of the experimental mass positions was suggested in 
Ref. [36], where a systematic lowering in mass of K°p peaks relative to the K + n was noted. In the 1 = 1 channel, 
there is no significant © ++ signal in the experimental data of the invariant K + p mass in the photo-induced reactions 
[5,9,37]. It is important that the widths of the J 77 = 3/2~ 0/=i should be about one order smaller than those of 
the J 77 = l/2 + and 3/2 + for the ©~j~. For the ++ search, it would be helpful to choose proper entrance and decay 
channels based on the further investigation of the production mechanism. In order to compare the present findings 
with the experimental data is more details, further experimental data with high resolution and high statistics are 
required. 

Finally, we would like to remind the readers that the absolute masses of the pentaquark in the present work 
are not predictions. We have an ambiguity of the zero-point energy of the string potential, which depends on the 
flux-tube topology in each of meson, three-quark baryon, pentaquark systems. In the present calculation of the 
pentaquark, we phenomenologically adjust it to reproduce the observed mass of the + . In order to predict absolute 
masses of unknown multiquarks with new flux-tube topologies, it is desirable to determine the zero-point energy more 
theoretically. 

The authors would like to thank to T. Kunihiro, Y. Akaishi and H. En'yo for valuable discussions. This work 
was supported by Japan Society for the Promotion of Science and Grants-in-Aid for Scientific Research of the Japan 
Ministry of Education, Science Sports, Culture, and Technology. 
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